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At a telephone exchange, calls appear before a simple trunk group of 
m lines in accordance with a recurrent process. If every line is busy, calls 
are delayed. The call holding times are mutually independent random 
variables with common exponential distribution. In this paper, methods are 
given for the determination of the distribution of the delay for a stationary 
process and various orders of service. Three orders of service are considered: 
(1) order of arrival, (2) random order, and (3) inverse order of arrival. 

I. INTRODUCTION 

In the theory of telephone traffic, the following process is of consider- 
able interest. In the time interval ^ I < <», calls appear before a 
simple trunk group with m lines at instants n , t 2 , ■ • ■ , t„ , • • • where 
the interarrival times r n+ i - rjn = 1,2, • • • ) are identically distrib- 
uted, mutually independent, positive random variables with distribu- 
tion function 

P{r n+1 - r„ :g x) = F(x) (n = 1,2, •■•)• (1) 

We say that the call input is a recurrent process. If an incoming call 
finds a free line, a connection is realized instantaneously. If every line 
is busy, the incoming call is delayed and waits for service as long as 
necessary (no defections). Denote by x« the holding time of the nth 
call. It is supposed that {x,.| is a sequence of identically distributed, 
mutually independent, positive random variables with distribution func- 
tion 

Pj x „ <k x\ = H(x) (n = 1,2, -..) (2) 



t An address presented on September 14, 1961, at the Troisieme Congrcs Inter- 
national de Tdletrafic, Paris. 



311 



312 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1902 

and independent of {t,,}. Wc shall consider only those systems of service 
which satisfy the requirements that there is no free line if there are calls 
waiting and that the same principle of service applies to every call (no 
priorities). Such a service system can be characterized by the symbol 
[F(x),H(x),m] provided that the order of service is specified. 

The ideal order of service, "order of arrival" or "first come — first 
served," is not always realizable, particularly at times of heavy traffic; 
therefore it is important to consider other orders of service also. One of 
these is "service in random order" which often describes the practical 
situation with high accuracy. In this case, waiting calls are chosen for 
service at random. Every call, independently of the others, and of its 
past delay, has the same probability of being chosen. Further, it is of 
great informative value to consider the extreme case, "inverse order of 
arrival," or "last come — first served." (At present we are not concerned 
with "priority systems" in which "last come — first served" service is 
the natural order, e.g., the last information to be received may be the 
most important in the process of the arrival of messages.) 

In what follows we shall consider the system [F(x),H(x),m] in the 
particular case when call holding times have the exponential distribu- 
tion 

fl - (T" 1 (x ^ 0) 
H(x) = - (3) 

(0 (x < 0) 

and the process is stationary. We shall give methods for finding delay 
distributions for the three service orders mentioned. 
We introduce the notation 

<p(s) = f e~"dF(x) (4) 

for the Laplace-Stieltjes transform of the distribution function of inter- 
arrival times, 

a = J xdF(x) (5) 

for the average interarrival time, W(x) for the delay distribution func- 
tion, i.e., W(x) is the probability that the delay is ^x, and 

q( 8 ) = f e' 8X dW(x) (6) 

for the Laplace-Stieltjes transform of W(x). 
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If 9t(s) ^ then denote by z = y(s) the root, with the smallest ab- 
solute value of the equation 

z = <p(s + mn{\ — z)). 

We have | y(s) | ^ 1 and for r = 1,2, • ■ • 

w^-'S^jif •"***" ^"-^ (7) 

where F n (x) denotes the rath iterated convolution of F(.r) with itself. 
Let co = t(0); then 

.-Z^fi-^ff.!.). (8) 

If /Hflfi ^ 1, then co = 1 while if man > 1, then co is real and < co < 1 . 

II. GENERAL THEORY 

A. K. Erlang 1 was the first to consider the process [F(x),H(x),m] in 
the particular case F(x) = 1 - e'^(x ^ 0), H(x) = 1 - e'" z (x ^ 0). 
The case of general F(x) has been treated earlier by D. G. Kendall," 
F. Pollaczek, 3 and the author. 4 

Denote by £(0 the number of calls in the system at the instant t; 
i.e., £(0 is the total number of calls either waiting or being served. 
Denote by x(0 the time difference between t and the arrival of the next 
call after t. Let £„ = £(t„ — 0), i.e., the nth call finds £„ calls in the sys- 
tem, and denote by 77,, the delay of the nth call. The initial state is given 
by ?(0) andx(0). 

The vector process |£(0>x(0; ^ t < 00) is a Markov process and 
has the same stochastic behavior for each order of service provided that 
there is no free line if there are calls waiting. In Ref. 4 it is proved that if 
man > 1> then there exists a unique stationary process. By choosing 
the suitable distribution for !£(0),x(0)| we arrive at the stationary proc- 
ess. For the stationary process, \k(t),x(t)\ has the same distribution 
for all /, and the distribution of £„ is independent of n. Let Pj£„ = k] = 
P k (k = 0, 1 , • • • ). As shown in Ref. 4 



1\ = 



E (-ir*n v r (k = 0, 1, ••• , m - 1) 

W (9) 

ylco '""' (k = m, m 4- 1, • • • ) 
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where 

U r = AC T t ( m ) ^V^'m""? , (10) 

s-ri i \-t( m ) fr* 1 -*<'"»- j r (ID 

A \l - « + M \jj Cy[l - *»(JM)][W(1 " -) - JV 

C '-U^m) ' (12) 

and co is defined by (8). 

Remark 1 — In the particular case when { r„} is a Poisson process of 
density X, i.e., F{x) = 1 - e~ Xz if x ^ 0, we have <p(s) = \/(X + s) 
and thus 



and 



A = 



CO = 


_X_ 
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Sjiw + »iw v «mJ 



Also for the stationary process the distribution function of ?/„ is inde- 
pendent of n but depends on the order of service. We shall use the nota- 
tion P{?j n ^ x\ = W(x) for all cases. It is to be noted that the expecta- 
tion E{ ??„} is independent of the order of service if the same principle of 
service applies to every call. We shall see later that in each case, the 
mean waiting time is given by 

f 'xdW(x) =-—^ r- 2 . (13) 

III. SERVICE IN ORDER OP ARRIVAL 

The following theorem has been proved earlier by D. G. Kendall, 
F. Pollaczek 3 and the author. 4 

Theorem 1 — The delay distribution function for order of arrival service 
is given by 

W (x) = 1 - t^— e-""^* (x ^ 0) (14) 

1 — CO 

where co is defined by (8) and A by (11). 
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Proof — We have for x ^ that 

1 — 1 « /*I / \fc 

/t=o fc=o ""o fc! 



lT(a:) = £ P* + Z P,„ +t f e-'""" ^^ m M d M (15) 



where {Pk\ is defined by (9). For, if an arriving call finds a free line, 
which has probability 

m— 1 o° a 

2^ Pk = 1 — t\j Pm+k = 1 — Z 

k=0 fc=o 1 — CO 

then its service starts without delay; if it finds every line busy and k 
(k = 0, 1, • • • ) calls waiting, which has probability P m+ k = Aco k , then 
its service starts at the (k + 1 )st departure after the arrival. Since the 
departures follow a Poisson process of density mn, under this condition, 
the probability that the delay ^ x is 

f* -m» u (mnu) k , 
/ e — j-: — mn au. 

•>o kl 

Thus (15) follows from these and (14) agrees with (15). 
It follows from (14) that the average delay is 

[ xdW(x) =—- (16) 

J mn{l — co) 1 

and the second moment of the delay is 

f x 2 dW(x) = ^ . (17) 

Jis (m,u) 2 (l — air 



mi 



IV. SERVICE IN RANDOM ORDER 



In the particular case of Poisson input this process has been investi- 
ated by S. D. Mellor, 5 E. Vaulot, 8 C. Palm, 7 ' 8 f F. Pollaczek, 9 J. Rior- 
dan, 10 R. I. Wilkinson, 11 and J. LeRoy. 12 Now we shall consider the case 
of recurrent input. 

Denote by Wj(x) (j = 0, 1, • • •) the probability that the delay of 
a call is ^.r, given that on its arrival all lines are busy and j other calls 
are waiting. For the stationary process the probability that an arriving 
call finds all m lines busy and j calls waiting is Pj+ m = Au) J . If there is 
a free line when a call arrives, which has probability 
"^} A 

Zft-i- 



,=o 1 — u 

then there is no delay. 

t Ref. 8 is an English version of the material in Ref . 7. 
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Thus we have 



that is, 



W(x) = £ P, + S P m+ jWj(x), (18) 

7=0 J=0 



W(x) = 1 - — ^— + 4Z lF y (aO«'. (19) 

1 — co y=o 



Let us introduce the notation 



Oi(e) = / e—'dWiix) (20) 

•'o 

and 

00 

•(m) = EoyCay (2i) 

7=0 

which is convergent if 9?(s) ^ and | 2 | < 1. 

Theorem 2 — The Laplace- Stieltjes transform of the delay distribution 
function for random service is given by 

0(a) = 1 - — h A$(s,to) (22) 

1 — O) 

where co = 7(0) is defined by (8), A by (11), and 

. , mfi f [' du 

5/05 = ; exD \ 1 

8 + m/*[l — 7(«)] \Jyi8)<p(s + mn(l — u)) — u 



+ WM / 7i 

J 7(«) (1 — 



[1 - <p(s + m/x(l — w))] 



(23) 



w)(s + w/i(l — u))[u — <p(s + m/i(l — w))] 

■ exp {/. rf. + mM(t -.))-. } d " 

where y(s) is the only root in z of the equation 

z = <p(s + m/x(l - 2)) (24) 

m f/ie wmY aVc/e | 2 | < 1. The explicit form of y(s) is given by (7) with 
r = 1. 

Proof — For j = 0, 1, • • • we can write that 



fc=0 7 + 1 J ft! 



(25) 
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where * denotes convolution. If an arriving call finds every line busy 
and j ' ( j ' = 0, 1 , • • • ) calls waiting, then the event that the delay is ^ x 
can occur in the following mutually exclusive ways: either in the sub- 
sequent interarrival interval k (k = 0, 1, • • • , j) services terminate 
and the service of the given call does not start during this time interval 
or, in the subsequent interarrival interval at least k -\- 1 (k = 0, 1, 
• • • , j) services end and the service of the given call starts at the termi- 
nation point of the (k + l)st service. 

Forming the Laplace-Stieltjes transform of (25) we get 

(i+l)0,(s) = ZO+ 1 - *)«*■!-*(«) (Y ( * w ^f^rfF(.r) 

fe-0 •'o A' . 

i - , ,* (26) 

+ Z / e-**-* ^L [1 - F(x)]mn dx 

whence 

00 00 

z E (./ + 1 )«/(«)*' = ¥>(s + mud - z)) T t jQ j (s)z i 

i=0 i—o 

muz [1 - y?(s + m/i(l - z))] 
"•" (1 - 2 ) (s + m M (l -2)) ' 

that is, 

[z - <p(s + mn(\ - z))} ^p^ + $(s,z) 

dz 

(27) 
/?(M [1 — <p($ + //nx(l — 2))] 



(1-2) (8 + «|i(l-a)) 

If m^« > I, then | p(s -f- m/z(l — z)) | 22 tp{mp.t) < 1 — e when | 2 | = 
1 — e and e is a sufficiently small positive number. Consequently by 
Rouch6's theorem it follows that 

2 = <p(s + mn(l - 2)) 

has one and only one root 2 = y(s) in the circle | z | < 1 — e, where 
e is a sufficiently small positive number. The explicit form (7) for 
[7(6') ) r can be obtained by Lagrange expansion. By definition &(s,z) is 
a regular function of 2 if | 2 | < 1 and SR(s) ^ 0. If we put 2 = y(s) in 
(27), then we get 

*(« >7 (*)) =— r — ?r — rvi- (28) 

s + w M [l - 7(*)1 

The solution of the differential equation (27) which satisfies (28) can 
be written in the form (23). Finally, (22) follows from (19). 
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Remark 2 — Let us introduce the notation 

A{s,z) =z- V {s+ m M (l - z)) (29) 

and 

D , v mn [1 - <p(s + mn{l - z))\ ,„ m 

Then (23) can be written in the following equivalent form 

$(s,z) = B(s,z) - exp<-/ — — r> — du. (31) 

The function $>(s,z) can also be expressed in the form of an infinite 
series as follows 

*(M) = Z^(s)[2-T(s)] J ' (32) 

j'=o 

which is convergent if | z — y(s) | is small enough. If 

A(s,z) = t,Aj(8)[z-y(8)) J (33) 

[note that A(s,t(s) ) = by definition of y(s)] and 

B(s,z) = T,Bj(8)[z-y(8)y (34) 

/=o 

then 4>;(s) (j = 0, 1, • • • ) can be obtained by the following recurrence 
formula 

£ feMs)Ay +1 -*(s) + *,■(«) = B,(«) (j = 0, 1, •••)• (35) 



<t=o 



This follows from (27). In particular by (35) we obtain 
«t>o(s) = B (s) , $i(s) = y-qrxTs)' 
5 2 (s) fii(s)^ 2 (s) 



* 2 (s) = 



1 + 2AM [1 + Ai(s)][l + 2A 2 (s)]" 

Formula (32) can conveniently be used to determine the moments of 
the distribution function W(x). The rth moment 

[ xdW{x) 

can be calculated by the aid of the derivatives 

(i +j £r). 



m 
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By using the relation 7 (s) = <p(s + m/n(l — y(s)) we can write that 
sip (wju(1 — w)) 



7(s) = w + 



[1 + mn<p'(mn(l — w))] 



_ s 2 ip"(mn (1 - a)) , 

2[1 + Wtod -c))] 3 "^ 



(36) 



Now in particular we have / x dW(x) = — 7- rr (37) 

Jo mull — o)j- 



and / x* dW(x) = —r - r — - w ^ . 

J (mn)-(l — co) 3 \_2 + mn<p'(mn(l — w))J 



(38) 



V. SERVICE IN INVERSE ORDER OF ARRIVAL 

The particular case when {t„\ is a Poisson process was investigated 
earlier by E. Vaulot, 13 J. Riordan, 14 and D. M. G. Wishart. 15 The case 
of recurrent input can be treated in a similar way. As noted by J. Rior- 
dan the problem can be reduced to rinding the distribution of the length 
of the busy period for the process of type [F(x),H(x),\] where F(x) is 
defined by ( 1 ) , 

and there is only one server. In this case denote by G(x) the probability 
that the length of the busy period is ^x. The busy period is defined as 
the time interval during which the server is continuously busy. Evi- 
dently G(x) is independent of the order of service, provided that the 
server is idle if and only if there is no waiting customer in the system. 

If mna > 1, then there is a unique stationary process, and for the 
stationary process W(x) is given by 

Theorem 4 — The delay distribution function for last-come, first-served 
service is 

W(x) = 1 - -A- + -A_ G(x) (40) 

I — CO 1 — 0} 

where co is defined by (8), A by (11), and for x ^ 

G(x) = mn £ e-" 1 (mfiX -" f [1 - F n (u)] du (41) 

n=l n\ Jq 

with F„(u) the nth iterated convolution of F(u) with itself. 

Proof — If a call arrives and finds a free line, which has probability 
1 — A(l — u)~ , then its service starts without delay; if on its arrival 
every line is busy, then we can remove all the calls waiting without 
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affecting the distribution of the delay of the call in question. The service 
of this call starts when the queue size decreases to m for the first time. 
The waiting time of this call evidently has the same distribution as the 
length of the busy period for the queueing process of type [F(x),H(x),l] 
with H(x) = 1 — e~ mi11 (x ^ 0). For, in both cases the arrivals have 
identical stochastic law and the departures follow a Poisson process of 
density mn. Thus we get (40). In Ref. 16 it is proved that 

■Jgfa)- ™ A1 ~ yis) ) , } (42) 

where y(s) is the root with smallest absolute value in z of the equation 

z = *(s + m M (l - *))• (43) 

y(s) is given by (7) with r = 1. By Lagrange expansion we find that 



/; 



ant \ m f i _i_ „V (-1)"(wm)" 

dG(x) = — +s2, —. 

s -f- nifj. h=i n\ 



d"- 1 ([<p{s + m»)] n 



ds"- 1 \ (s 4- ^m) ! 



(44) 



whence (41) follows by inversion. 
By using the expansion (36) we get from (42) that 



and 



f 

•'u 



f x dW(x) = —-^ ra (45) 



24 
X 2 dW(x) = 7 ttt-j wT~r 77 71 TVi- ^^ 

(mn)-(l — w) 3 [l + mw> r (mp(i — w))l 
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